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If^ ' In this paper we investigate a parametric resonance phenomenon of a 

o ■ 

Kaluza-Klein mode in a D-dimensional generalized Kaluza-Klein theory. As 
I the origin of the parametric resonance we consider a small oscillation of a 



scale of the compactification around a today's value of it. To make our ar- 



^j-jl guments definite and for simplicity we consider two classes of models of the 



compactification: those hy Sd {d = D — 4) and those by Sd^ x Sd2 {di > d' 



2, 



?H ' di + d2 = D — 4). For these models we show that parametric resonance can 

occur for the Kaluza-Klein mode. After that, we give formulas of a creation 

rate and a number of created quanta of the Kaluza-Klein mode due to the 

parametric resonance, taking into account the first and the second resonance 

band. By using the formulas we calculate those quantities for each model of 

the compactification. Finally we give conditions for the parametric resonance 

to be efficient and discuss cosmological implications. 
PACS numbers: 04.50.-hh, 98.80.Cq 
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I. INTRODUCTION 



Many unified theories require spacetime dimension higher than 4. We call them general- 
ized Kaluza-Klein theories. For example the superstring theory predicts spacetime dimension 
-D = 10 [ip, the M-theory D = 11 0. On the contrary we can see only a 4-dimensional part 
of the spacetime by experiments. Hence the extra {D — 4)-dimensional part is thought to 
be compactified so that it cannot be seen by us. One may be able to give a 4-dimensional 
theory describing our 4-dimensional universe from a generalized Kaluza-Klein theory by the 
compactification of the extra dimensions. However, there are a plenty number of ways of 
compactification and there are the same number of the corresponding 4-dimensional theo- 
ries for a D-dimensional theory. Thus we need criterion to judge which compactification 
must be selected. For example, in the case of the superstring theory it is well known that a 
'good' compactification must be realized by a Calabi-Yau manifold in order for the = 1 
supersymmetry to remain in the corresponding 4-dimensional theory In this paper we 
give a suggestive argument toward a cosmological criterion for compactification. 

Kolb and Slansky examined a model of the compactification by 5*1 in the 5-dimensional 
theory with a metric field described by an Einstein-Hilbert action and a massless scalar 
field described by a Klein-Gordon action. They investigated a cosmological evolution of 
energy density of a mode of the scalar field with a non-zero momentum in the direction of 
the compactified extra space Si assuming that there is no entropy production. The mode 
is usually called a Kaluza-Klein mode and the quantum of the mode is called a pyrgon. 
They showed that if a scale of the compactification is of order of the 4-dimensional Planck 
length and energy density of pyrgons is comparable with energy density of radiation at an 
early epoch of the universe, then at the recombination epoch the former energy density 
terribly exceeds the critical density of the universe. The physical reason of their result is 
essentially the momentum conservation in the direction of the compactification: a pyrgon 
with a positive momentum in the direction cannot decay without meeting with another 
pyrgon with a negative momentum. For this reason energy density pkk of the Kaluza- 
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Klein mode decreases as a~^, where a is a scale factor of the universe, while energy density 
Prad of radiation decreases as a""^. As a result the ratio Pkk/ Prad evolves as a. Since the 
momentum conservation in the direction of the compactification does always hold, their 
result is expected to hold in a wider class of models of the compactification. Thus, if 
entropy production is negligible, energy density of pyrgons must not be comparable with 
energy density of radiation at any early epoch of the universe in a wider class of models 
of the compactification. On the other hand, if pyrgons are suitably created at an early 
epoch and they are diluted by entropy production after their creation (eg. reheating after 
inflation) then the pyrgons may become a dark matter at today. Anyway, we have to analyze 
whether pyrgons can be created in early universe or not in order to obtain a knowledge of 
the compactification from a cosmological point of view. In this paper, motivated by these 
considerations, we investigate a catastrophic creation of quanta of a Kaluza-Klein mode in 
a D-dimensional generalized Kaluza-Klein theory. 

As the origin of the catastrophic creation of quanta of the Kaluza-Klein mode we consider 
a small oscillation of a scale of the compactification at an early epoch of the universe. In 
general it is expected that frequency of the oscillation is of order of inverse of the scale 
of the compactification and mass of the Kaluza-Klein mode is of the same order. Hence 
the mass of the Kaluza-Klein mode oscillates with a frequency of order of the mass itself. 
Such an oscillation of the mass causes a catastrophic creation of quanta when frequency and 
amplitude of the oscillation are in a band in their parameter space. This sort of catastrophic 
creation of quanta is called a parametric resonance phenomenon [Q. In this paper we apply a 
theory of parametric resonance [| to a Kaluza-Klein mode and judge whether the catastrophic 
creation occurs efficiently or not. 

In Sec. we calculate 4-dimensional actions for two classes of simple models of the 



^ The theory of parametric resonance was applied to reheating (or preheating) process after 
inflation by many authors ||5|,^. 
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compactification. In Sec. |T| a Kaluza-Klein mode is investigated quantum mechanically 
on a classical background of a 4-dimensional Friedmann universe with a small oscillation of 
the scale of the compactification. We show that parametric resonance of the Kaluza-Klein 
mode can occur. After that we intend to judge whether the resonance is efficient to produce 
sufficient quanta of the Kaluza-Klein mode. Sec. is devoted to summarize this paper. 



II. 4-DIMENSIONAL ACTION 

Our starting point in this paper is a D-dimensional Einstein-Hilbert action with a cos- 
mological constant A: 

lEH = ^j d'^xV^iR - 2A), (2.1) 

where k is a positive constant, g and R are a determinant and a Ricci scalar of a D- 
dimensional metric tensor Qmn Q- In addition to the gravitational field described by the 
Einstein-Hilbert action, we consider other massless fields, too. For the present, we include 
a real massless scalar field (p described by a /^-dimensional Klein-Gordon action 

Ikg = -\j d^'x^g^'^'dM^dN^. (2.2) 

As stated in the previous section, some candidates of the unified theory require D > A while 
the spacetime dimension we can see by experiments is 4. Hence we have to compactify the 
extra {D — 4)-dimensional space so that it can not be seen by experiments. In this paper to 
make arguments definite we consider two classes of simple models of the compactification. 
Unfortunately, there are a great number of other models of the compactification and there is 



^ Hereafter we take the sign convention ( — !- + •••+) of the metric tensor. 

^ Here we mention that the low energy action of the superstring theory and the M-theory can 
be reduced to a form which includes the Einstein-Hilbert action and the massless Klein-Gordon 
action. 
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no conclusive first principle to determine which model should be taken. However, qualitative 
properties shown in this paper for simple models are expected to be common in a wider class 
of models of the compactification. 



A. Compactification by 

First we consider a compactification by a d-dimensional sphere Sd {d — D — A). In this 
case the metric tensor gMN can be written as 

guNdx^dx^ = g^.^dx^dx" + b^Q\fdx'dx^, (2.3) 

where is a 4-dimensional metric depending only on the 4-dimensional coordinates x'* 
(// = 0, 1, 2, 3), 6 is a scale of the compactification (a radius of the d-dimensional sphere) 
depending only on the 4-dimensional coordinates, and ^Ifdx^dx^ is a line element of a unit 
d-sphere. The scalar field 4> is expanded as follows: 

l,m 

where the constant bo is a today's value of b. In the expression, F^^-* {I — 1, 2, • • •; and m 
denotes a set of d — 1 numbers that are need in order for a set of all YJ^-* to be a complete 
set of functions on the (i-sphere) are real harmonics on the d-sphere satisfying 

-^d, (VliWQW'^-S^.y +l{l + d- 1)1^^) = 0, (2.5) 

J d'^xVmYl-yit' = ^ll'^mm'. (2.6) 

and (f)im{x^) is a real function depending only on the 4-dimensional coordinates x^. The 
Einstein-Hilbert action and the Klein-Gordon action in this case are 

d 



Ieh = ^ /"^'^V^d;) [R + d{d-i)r''{d^\^b){dAnb) + d{d-i)b-^ -2K 



(2.7) 
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where k is a positive constant defined by — / (2^6^ tt) and ^ is a Ricci scalar of the 
4-dimensional metric tensor Since we shall consider the case when the scale h of the 
compactification is time dependent, the factor (b/bo)^ before the Ricci scalar R and before 
the kinetic term g^"d^ji(t)ijndv(t>im of the field is not a constant. Hence those expressions 
of the Einstein-Hilbert action and the Klein-Gordon action are quite different from the 
usual form of the 4-dimensional Einstein-Hilbert action and the 4-dimensional Klein-Gordon 
action. 

In this paper we shall examine excitations of the 0/^ fields due to a time dependence 
of the compactification scale b. For this purpose we have to analyze a quantum dynamics 
of the fields. However, as seen in the previous paragraph, the action of the (f)im fields 
looks quite different from the usual 4-dimensional Klein-Gordon action and is not a useful 
form for the analysis of the quantum dynamics. Thus we perform the following conformal 
transformation in order to obtain a useful expression. 

g„v = g„v (2.8) 

By this conformal transformation we obtain 

I EH = j d^xy/^ —R - -g'^'d^ad^a - Uo{a) , 

Ikg = E / d'^V^ [g'^^d.^PiMim + M2(a)0L] , (2-9) 
where is a Ricci scalar of the 4-dimensional metric tensor 51^^,, the field a is defined by 



..W^, (.10, 



the mass Mi{a) of the 4-dimensional scalar field (j)im is given by 



Ml{a) = ^) ^-{d+2)a/a,^ ^2.11) 

^0 



and the potential Uq^u) of cr is 



Uo{a) = Ae-'i-M _ ^((^ J) ^-id+2)a/.o_ (2.12) 

Here we mention that a mode with non-zero / has a mass of order I/60 and is called a 
Kaluza-Klein mode. 

The expression ( p^.9| ) can be understood as the 4-dimensional Einstein-Hilbert action plus 
the 4-dimensional Klein-Gordon action with a potential. However, it can be easily confirmed 
that UQ{(y) has no local minimum. It means that there is no stable compactification by 
Sd if the model is not modified any more. Therefore we have to take into account some 
mechanism in order to stabilize the potential of the a field. In this paper as the mechanism 
we consider a 1-loop effective action (an action for so-called Casimir effects) contributed by 
all 4-dimensional matter fields. The 1-loop effective action is of the following form 0: 

hioop = - j d^x^/^Viioop{,CF), (2.13) 

where Vizoop is an effective potential, which in turn is a sum over all contributing fields 0. 

Viioop{(y) = Y.c-y{M,), (2.14) 

i 

where i distinguishes 4-dimensional fields contributing to the 1-loop effective action, is 
a mass of the field specified by and Cj is a numerical factor depending on what spin the 
field i has 0. In the expression the function V{M) is given by 

ViM) = -- lim(47r)-"/2r(_n/2)M". (2.15) 

2 n^4 

Since squared masses of 4-dimensional fields deduced from the massless D-dimensional the- 
ory is proportional to e~('^+^)'^/'^° as the squared mass ( p.ll|) of the (pim fields, the 1-loop 
effective action is of the following form in general. 



^ Strictly speaking, we have to restrict the dimension D to be odd since otherwise the conformal 
anomaly arises and the arguments become complicated. 
^ For example, Cj = 1 for spin 0, Cj = —4 for spin 1/2 0. 



V^ 



Hoop 



(J 



^g-2(d+2)(T/o 



(2.16) 



where A is a constant. We could calculate the constant A in principle if we knew details 
of the D-dimensional theory. However, we do not know it. Thus we shall determine the 
constant A phenomenologically. We require that the total potential f/o(cr) + Vnoopicr) of the 
field a has an extremum at a = (or at 6 = h^, where 6o is the today's value of h) and the 
extremum (or the 4-dimensional cosmological constant) is zero. Thus the constant A and 
the renormalized value of the dimensional cosmological constant A can be expressed in 
terms of h^, and the total potential of a is 



a 



_d + 2 

where a is a constant defined by 



g-2(d+2)(T/o-o _j_ g-dcr/o-o 



d + 4 
d + 2 



~id+2)a-/ao 



(2.17) 



a 



d{d- l)(d + 2) 



2{d + A)K%l ' 

This form of the total potential of the field a is used in many literature 
derivative at cr = is given by 

2{d-l) 



1? 



(2.18) 
Its second 

(2.19) 



and is positive for d > 2. Thus, for d > 2, the potential Ui has a local minimum at o" = 0. 
Finally the 4-dimensional action for the fields gfj,,^, a and (j)im in this model is 



d'^x 
1 



l,m 



EM' 
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(2.20) 



B. Compactification by Sd^ x 

Next we consider a compactification by a direct product 5*^^ x 5*^2 of a rfi-dimensional 
sphere and a (i2-dimensional sphere {di + d2 = D — 4) . We assume that di > d2 > 1. In this 
case the metric tensor qmn can be written as 
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gMNdx^dx^ = g^^dx^'dx" + bjn\f^dx'dx^ + bln^^.^^dx^dx'^, 



(2.21) 



where ^^i, is a 4- dimensional metric depending only on the 4-dimensional coordinates x'^ 
{fi = 0,1,2,3), bi and 62 are radii of the (ii-dimensional sphere and the (i2-dimensional 
sphere depending only on the 4-dimensional coordinates, and Q^j^^dx^dx^ and Q^^^^'dx^dx'^ 
are line elements of a unit rfi-sphere and a unit (i2-sphere, respectively. The scalar field is 
expanded as follows: 



,-di/2,-d2/2 
"10 ^20 



0/l/2mim2 i^'^)^hm2 (^*)^2m2 (^^) 

Il,l2,mi,m2 



[2.22) 



where 610 and 620 are today's values of bi and 62, ^ifm^ and Yi^^^ are real harmonics on 
the (ii-sphere and the (i2-sphere, respectively, and 4>iii2mim2{x^) is a real function depending 
only on the 4-dimensional coordinates x^. The Einstein-Hilbert action and the Klein-Gordon 
action in this case are 



d'^x^—g 
1 



-^7? - ]-g^'''d^a+d^a+ - ^g^'^d^a^d^a^ - Uoi^i, ^2) 



IkG = -7; d^Xy^ g^''d^(j)i^i^rnirn2du(Phl2m^m2 + M^,l2i(^l,(^2)4>ll2 



m\m2 



(2.23) 



I\,l2,m\,m2 ' 

where we have performed the conformal transformation Q 



9^u 



bi\ / 62 1 



(2.24) 



,^10/ \b20 ) 

In the expressions, k is a positive constant defined by = / {2'^^^'^'^biQb'^Ti'^) , R is a Ricci 
scalar of the 4-dimensional metric tensor g^^,. The 4-dimensional fields ai^2 are defined by 

^ ^1,2 \ 



o"i,2 = crio,20 In 



^10,20, 



C"10,20 



l di,2{di,2 + 2) 
2h? 



(2.25) 



and (j± are linear combinations of ai 2: 



^ See the comments before (|2.8| 
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^± = y^(^i±^2), (2.26) 

where a = \J (^di+%(d^ +2)- "^^^ mass Mi^^i^^cri, 0-2) of the 4- dimensional scalar field 4>i^i2mim2 is 
given by 



'^lO ^^20 



(2.27) 

and the potential t/o(cri,a2) is 



Uo{(ri, 0-2) = e" 



-dlCTi /(TiO — d2O'2/o'20 



A_ _ rfl(rfl-l) 2.,Mo _ 4(4-1) 2.2M0' 
^2 2/€262 2/t262 



(2.28) 



Here we mention that a mode with non-zero li or non-zero ^2 has a mass of order I/610 or 
1/620 and is called a Kaluza-Klein mode. 

As in the case of the compactification by Sd, it can be easily confirmed that f/o(o"i)Cr2) 
has no local minimum, and that there is no stable compactification by 5*^^ x 5*^2 if the model 
is not modified any more. In order to stabilize the potential we consider a 1-loop effective 
action (an action for so-called Casimir effects) contributed by all 4-dimensional matter fields. 
The 1-loop effective action is of the following form |^: 

hioop = - J d'^x^/^Viioop{,(yi,(r2), (2.29) 
where Vuoop is an effective potential, which in turn is a sum over all contributing fields []. 

ViiooM^^2) = Ec.^(^.), (2.30) 

i 

where i distinguishes 4-dimensional fields contributing to the 1-loop effective action. Mi is 
a mass of the field specified by and q is a numerical factor depending on what spin the 
field i has 0. In the expression the function V{M) is given by ( |2.15| ). Since, in the present 



^ See the footnote soon after ( p.l3| ). 
See the footnote soon after (|]l|). 
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model, squared masses of 4-dimensional fields deduced from the D-dimensional massless 
theory is a sum of a term proportional to e~'^'^i+^^°"i/°""'~°'2'^2/o'2o g^j^^ term proportional to 
g-diai/,7io-(d2+2)a2/a2o gg ti^g squared mass (^27^ of the (pim fields, the 1-loop effective action 
is of the following form in general. 

(2.31) 

where An, A12 and A22 are constants. We require that the total potential Uo{ai,a2) + 
Viioop{cri, o'2) of the fields di and a2 has an extremum at ai = (72 = (or at 61 = 610, &2 = 
where 610 and 620 are the today's values of bi and 62? respectively) and the extremum (or 
the 4-dimensional cosmological constant) is zero. Thus An, A22 and A can be expressed in 
terms of 610, &20 and the unknown constant A12, and the total potential of cxi and a2 is 

^^1(0-1,0-2) = aiUi^\ai,a2) + a2f^fVi>o-2) 

j^_^^Q — 2dicri/(7io—2d2Cr2/cr20 ^g — 2o-i/ctio g— 2CT2/f20 ^ ^2 32) 

(1 2) 

where the constants ai^2 and the functions f/j ' {0-1,0-2) are defined by 

dijdi - l){di + d2 + 2) 



a2 



2{di + d2 + ^)ti^bjo 
d2{d2-l){di + d2 + 2] 
2{di + d2 + 4:)K% 



'20 



/-/-(l)/^ ^ ^ — d2 + 4: ^-2(di+2)cri/aio-2d2cr2/cr20 _|_ „-dio-i /(Tio- (i2 0-2/0-20 

^^^'^^^"2(rfi + d2 + 2)' 



(il + (i2 + 4 _(di+2)oi/oio-d202/o20 



-e 



di + ^2 + 2 

-2di 01 /oio -2(c(2+2)o2 /0-20 



^2 



2(^1 + c/2 + 2) 

7-^(2)/ n- — -2dioi/oio-2(ii2+2)o-2/o-20 i p -oil 01/010 -0(202/020 

2((ii + d2 + 2) 



dl + d2+4: _-diOi/oio-(d2+2)o2/o20 



-e 



+ ^2 + 2 

^ -2(di+2)oi/oio-2d202/o2o oq\ 

2(rfi + d2 + 2)" • ^^■'^'^^ 

Finally the 4-dimensional action for the fields f?^,^, Oi, 02 and </>iii2mim2 this model is 
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1 



mim2 ■ (2.34) 

Unfortunately, the potential f/i(c'"i, (T2) depends on the unknown constant A12. In order 
to make our arguments independent of the value of A12 we take the following ansatz: 



hi _ &2 
^10 ^20 



(2.35) 



Physically, this condition means that the shape of the extra [D — 4)-dimensional space 
Sd^ X Sd2 is fixed and only its volume changes. In this case let us introduce a field a by 



(t/cTo = (Ji/cTiQ = (J2/a20, 



'{di + d2){di + d2 + 2) 
2^2 



(2.36) 



The action ( p. 341) is simplified as 



1=1 d^xy/—g 
1 



h,l2,mi,m2 

where the mass Mi^^i^^a) of the Kaluza-Klein mode 0iii2mim2 is given by 

'h{h + di-l) 12(12 + d2- I] 



m2 



(2.37) 



/)2 
"20 



-{di+d2+2)cj/ao 



(2.38) 



and the potential Ui{a) of the a field is 



f/i(cr) = («! + 02) 



C?l + ^2 + 2 



_g-2(di+d2+2)o-/cro _|_ g-(dl+rf2)cr/o-o 



C?i + ^2 + 2 



(2.39) 



Note that the potential Ui{a) is independent of the unknown constant yli2.The second 
derivative of UAa) at a = is 



di{di - 1) d2{d2 - 1) 



(2.40) 



di + d2 

Hence, if di > 2, the extremum of Ui{a) at a = is a local minimum and the compactifica- 
tion is stable at least locally. 
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III. PARAMETRIC RESONANCE 



In the previous section we have shown that the compactification by Sd {d > 2) is stable, 
and the compactification by Sd-^ x {di > 2) is also stable at least in the sub-configuration- 
space defined by ( |2.35| ). Thus, we consider the Sd compactification with d > 2 and the 
Sdi X Sd2 compactification with di > 2. In this section we judge whether a catastrophic 
creation of quanta of the Kaluza-Klein mode (pim (or (t)i^i2mim2) can occur due to an oscillation 
of the a field around the local minimum a = 0. We treat the metric tensor g^^^ and the field 
a classically, and the fields (pim (or (t)i^i2m^m2) quantum mechanically. 



A. The hamiltonians for and 4)i^i2 



in\m2 



We assume that the 4-dimensional metric g^^ is of the Friedmann universe and that the 
field (J is homogeneous on this background. The following equation of motion for the field 
a is derived from the 4-dimensional action ( 2.20 ) or ( 2.37| ): 



a + ?>Ha + U[{a) = {), (3.1) 

where H = a/a (a is a scale factor of the universe) and the overdot denotes derivation with 
respect to the cosmological time t. Hereafter we assume that 

if<cu, (3.2) 

where uo is defined by 

^ = \\lum- (3.3) 

Hence for a small deviation from the local minimum a = (|cr/cro| ^ V-^)) ^ solution ij{t) 
of the equation ( pj.lD can be written as 

^ = 5-(t)cos[2c<;(t-to)], (3.4) 

where the coefficient a{t) (|cr(t)| ^ 1/-D) is a slowly varying function of t: 
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cr 

a 



(3.5) 



and to is an arbitrary constant. In order to determine a time dependence of the function 
5"(t), we introduce energy density p„ of the a field by 



From (|3.1| ) the following equation is easily derived: 



(3.6) 



= -3i7dr2. 



(3.7) 



By averaging this equation from t—TT/{2uj) to t+TT/{2uj), we obtain the following approximate 
equation. 



Since in our situation the value of is approximately given by 

1 



(3.8) 



P<x ^ |(2a;)2(aoa)2, 



the time dependence of a is derived from the averaged equation (|3.8| ) as 



(3.9) 



cr(t) oc a"^/2. 



(3.10) 



In the remaining of this section we investigate a quantum dynamics of the Kaluza-Klein 
mode (pim (for the Sd case) and (t)i^i2mim2 (for the 5*^^ x 5*^2 case) on the classical background 



of the Friedmann universe and the field a given by (|3.4| ) and (|3.1CI|) . For this purpose we 
first go to the momentum representation by 



0/^(x,t) = a 



'3/2 



(27r)3 



it) cos(k-x)+0S2'^(t) sin(k-x 



(3.11) 



iTmima (^) cos(k ■ x) + ^J.^f^.^^ (t) sin(k ■ x) , (3.12) 



14 



where (l)il^^^(t) and 0[j\'fmim2(^) ^^^^ functions of the cosmological time t, the vector x 
denotes the comoving coordinates f\. Next we obtain the Hamiltonian corresponding to the 
time t from the action (|2.2CI|) and ( p.37|) . The result is summarized as follows: 



Lm 



(2)k 
lm 



(3.13) 



for the Sri case, and 



(2)k 

i2mim2 



(3.14) 



h,h,mi,m2 ' 

for the S'di x 3^2 case. In the expressions the partial hamiltonians H^^^^ and h\^I^J^^^_^ are 
of the following form. 



p2 + ^]2(t)g2 



(3.15) 



where Q denotes (f)^^^^ for Z/;^^'^^'' and ^{^^'J,^^ ^S2m2. and P denotes its conjugate 
momentum. The function VL{t) is a positive function of t given by 



^^{1) =u? A + Be 



-e cos 2ujt 



(3.16) 



where the parameters A, B, e and are given by 

, 262 (k/a)' , , 



P 



d- 1 
e = (c/ + 2)a- 



c/- 1 

^ l{l + d-l), 



0( 



1 rf-1 



(3.17) 



^ To obtain the representation, for concreteness, we assume that the Friedmann universe is flat 
one. If we assume the open or closed Friedmann universe then the corresponding representation 



is changed. However, the resulting partial hamiltonians are of the same form as (3.15), provided 
that is replaced by the corresponding eigen value. 
We can set to = without loss of generality. 
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for the Sd case 0, and 



A 



2(cii + rf2)&?o(k/«) 



+ AA, 



B 



rfl((il-l)+rf2(rf2-l)(&10/M' 
2{di + d2) [h{ll +di-l)+ kik + d2- l)(6l0/&20)'] 



di{di-l) + d2{d2-l){bio/b2oy 
e = {di + d2 + 2)cr, 



\ 2{di + d2 



di{di ^ 1) _^ d2{d2 - 1) 



62 



10 



6? 



20 



(3.18) 



for the S'dj x 5*^2 case 0, respectively. In both expressions of A, the term AA denotes 



AA = -Lu-^ f^H' + 



(3.19) 



In the remaining of this paper we consider the case when the condition (|3.2|) and the 
following condition hold. 



(3.20) 



Hence we can neglect the term AA in A. In general the parameters A and e appearing in 
the function Vt have time dependence even when we neglect the term AA. However, since 
A oc and e oc a"^/^, it is shown from ( p.2| ) that their change is extremely slow: 



2H < UJ, 



-H « 



(3.21) 



If / = 0, then B = Q and the parametric resonance does not arise. Thus, hereafter we consider 
a mode with / > 1. 

^2 If li = I2 = 0, then B = and the parametric resonance does not arise. Thus, hereafter we 
consider a mode with /i > 1 or Z2 ^ 1- 
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B. The resonance band and the number of created quanta 



A system with a hamiltonian Hke ( 3.15| ), when it is canonicaUy quantized, has the prop- 
erty that a catastrophic creation of quanta (parametric resonance) can occur if the param- 
eters A, B and e is in a band in their parameter space. Our task now is to derive the band 
and a number of created quanta of the Kaluza-Klein mode (pim or (t)i^i2mrm2- First we expand 
both in power of e and in Fourier series. For this purpose we use the identity 



-e cos 2z 



/o(e) + 2£(-irJ„(e) cos2nz, 



n=l 



where /„(x) is the modified Bessel function and is expanded as 
The obtained expression is 



oo oo 



s=l n=— oo 

where z = cut. In the expression the constants h and g!^^ are given by 



(3.22) 



(3.23) 



(3.24) 



(n+2s) 
9±2n 



h = A + BIo{e), 

(_l)n+l5 



2"+2^s!(n + s)!' 

where tt, = 1, 2, ■ ■ ■; s = 0, 1, ■ ■ ■, and other g^^ are all zero. For example, 



B 



9^1 = -^, = (n^±4). 



B 



(3.25) 



(3.26) 



Next let us introduce a parameter A by 



P 



(3.27) 



where p and q are mutually prime positive integers. In appendix ^ and it is shown that 
the parametric resonance does occur if the parameter A is in the following band: 
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Ao-|5| < A< Ao+|(5|, (3.28) 
where Aq is a combination of the coefficients g^'' being of order 0(e), and S is defined by 

(3.29) 



92p/q 



E yn y2p/q-n (2) 



Thus the existence of non-trivial interval of the resonance band requires that 5 7^ 0, which 
in turn requires that p/g is 1 or 2. Physically the first resonance corresponds to a decay of 
a quantum of the field a with energy 2uj into two quanta of the Kaluza-Klein mode with 
each energy tu, and the second resonance corresponds to a decay of a quantum of a into one 
quantum of the Kaluza-Klein mode with energy 2uj. The parameters Aq and 5 in our case 
are given as follows: 

= - 32 ^' 

5 = I (3.30) 

for the p/q = 1 resonance; and 

= 24 ^' 

S = e (3.31) 

for the p/q = 2 resonance. Next the Floquet index which represents a 'strength' of the 
resonance, is given by the formula ( |A2(j| ) and can be expressed as 

= J^v/r^(l + 0(e)), (3.32) 

where x is defined by x = (A — Aq)/6. 

As shown in appendix by using the form of a number of excited quanta is given 
by the following formula: 

N{t) ~ sinh^ (^J fx+dz^ . (3.33) 
We now estimate the integral in this formula, li B < {p/qY then 



lim X = ±00, lim x = =Foo, 



(3.34) 



and thus we can use the following approximation: 



\e6\ 


dz 




2p/q 


dx 


I 

A=Ao 



J Vl — x'^dx 



|e5|7r 



dz 



dx 



(3.35) 



A=Ao 



The quantity |(i^/(ix|A=Ao ^e calculated from ( |3.3(J| ) or ( |3.31| ) as follows Q: 



dz 



dx 



A=An 



5 



4(5-1) 



Huj 



-1 



l + 0(e^ 



(3.36) 



for the p/q = 1 resonance with B < 1; and 



dz 



dx 



A=Ao 



B{B - 2) 



32(5-4) 



-{l + Oie 



(3.37) 



for the p/q = 2 resonance with B < 4. If 5 = 1 then x for the p/q = 1 resonance satisfies 

lim x = 0, (3.38) 

a— >{X> 

and thus the integral is estimated as 

fi^dz — ^ J ^dt^ (3.39) 

where t* is the time when the mode enters the p/q = 1 resonance band. If 5 = 4 then x for 
the p/q = 2 resonance satisfies 



lim X = 00, 



(3.40) 



and thus the integral / fj,^dz is given by (|3.35| ) at an order estimate, where \dz/dx\/\=Ao is 
given by 



dz 



dx 



3 1 



A=Ao 



2 Hco- 



(3.41) 



Therefore the number of created quanta by the p/q = 1 resonance is given by 



13 



To obtain the expressions we have used the fact that A oc a ^ and that e oc a 



-3/2 
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N ~ sinh 



25(5-1) ' Huj-^ 



(3.42) 



for B < 1; and 



N ~ sinh^ 



- / edt 

4 



(3.43) 



for i? = 1. Note that the integral in the expression ( |3.43 ) can be performed when a time 
dependence of the scale factor a is specified. The number of created quanta by the p/q = 2 
resonance is 

'B\B-2fTT 



N ~ sinh 



212(5-4) Hlu-^ 



(3.44) 



for 5 < 4; and 



~ sinh 



0(1) 



(3.45) 



for B = A. 

Finally we mention that in order for the parametric resonance to be efficient to produce 
the quanta of the Kaluza-Klein mode the following condition is necessary and sufficient: 



> 1, T>3H, 



(3.46) 



where F is a creation rate of the quanta at the center of the resonance band given by 



T = 2u fi+ 



(3.47) 



Using the formula ( p.32| ), F can be written explicitly as follows: 



V = -eu 



(3.48) 



for B < 1\ and 



B\B-2\ 
32 



(3.49) 



for 5 < 4. 
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C. The Sd case 



We now apply the results obtained in the previous subsection to each model of the 
compactification. First, in this subsection we consider the Sd case with d > 2. In this case 
it can be seen from ( |3.17| ) that B > 1 for / > 1 Hence the resonance of p/q = 1 does not 
occur. On the other hand, the resonance of p/q = 2 can occur for I = 1 since B\i=i = 4 for 
d = 2 and 2 < B\i=i < 4 for d > 3. For / > 2, the resonance of p/q = 2 does not arise since 
B> 4. 

For d = 2 and / = 1, the number of created quanta and the creation rate F are given 
by ( p.45| ) and (|3.49| ) with B = A. Hence the parametric resonance is efficient to produce the 



quanta of the corresponding Kaluza-Klein mode (pim when 

^>0(1)-^. (3.50) 
Po H 

where is energy density of the oscillation of a, which is approximately given by (|3.9|) , and 
Po is so called critical density of the universe defined by 

Po = (3.51) 



For > 3 and / = 1, and F are given by (^A§ and (ICTj ) with B = 2d/{d - 1) 



Hence the parametric resonance is efficient to produce the quanta of the corresponding 
Kaluza-Klein mode (pim when 

^>0(l).(|)^^ (3.52) 
Note that, because of the assumption ( p.2|) , the condition ( p. 501 ) or ( p.52|) requires 



^ > 1. (3.53) 
Po 

Thus in our situation {H <C uj and |cr/(To ^ V-^) can conclude that the parametric 
resonance do not overproduce the quanta of the Kaluza-Klein mode in this model. 



See the footnote soon after (|3l^) . 
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D. The X case 



Following the case, we consider the S^^ x 3^2 case. In this model the parameter B is 
given by ( |3.18| ). 

For li and I2 satisfying B < 1, the number N of created quanta of the corresponding 
mode and the creation rate F are given by ( |3.42| ) and ( |3.48| ), respectively. Hence, if i? < 1 
and l-B — 1| = 0(1), the parametric resonance is efficient to produce the quanta of the 
Kaluza-Klein mode (pi^i^m^m-z only when 

^ > 0(1) ■ ^. (3.54) 
Po H 

If i? < 1 and |i? — 1| = 0(e) then the condition for efficient production becomes 

^>0(1). (3.55) 
Po 

For li and I2 satisfying B = 1, N and F are given by (|3.43|) and (|3.48| ) with B = 1. In 
order to estimate the integral in the expression of we have to specify a time dependence 
of the scale factor a. If a oc c^* then 



I" 



~ , 2 e 
edt = - ■ — 
3 H 



(3.56) 

i=i. 



If a OC r (n > 2/3) then 



°° , 2n e 

eat = ■ — 

3n - 2 H 



(3.57) 



If a oc r (n < 2/3) then 



00 



edt = 00. (3.58) 



Anyway, the parametric resonance is efficient to produce the quanta of the corresponding 
Kaluza-Klein mode (pi^i^mimz when 

^>0(1). (3.59) 
Po 
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For /i and h satisfying 1<B<2ot2<B<4:, N and T are given by ( p.44|) and ( p. 491) 



respectively. Hence in this case, if |-B — 4| = 0(1), the parametric resonance is efficient to 
produce the quanta of the corresponding Kaluza-Klein mode 4>i^i2mim2 only when 

E.>0(l).{^f. (3.60) 



For li and I2 satisfying B = 4, N and F are given by ( p.45D and (|3.49| ) with B = 4. 



Hence the parametric resonance is efficient to produce the quanta of the corresponding 
Kaluza-Klein mode 0iii2mim2 when 

^>0(1)-^. (3.61) 

Finally for li and I2 satisfying B = 2 01 B > 4, N and F are zero in our treatment. Thus, 
in this case, there is no regime when the the parametric resonance is efficient to produce the 
quanta of the corresponding Kaluza-Klein mode (pi^i2mim2- 

After all we can conclude that if there is a mode such that -B < 1 and \B — 1| = 0(e) 
then the quanta of the Kaluza-Klein mode can be overproduced by the parametric resonance. 
Otherwise, at least in our situation {H ^ u and |a/cro ^ 1/-D) the parametric resonance is 
so mild that the overproduction does not occur. 



IV. SUMMARY AND DISCUSSION 

In this paper we have investigated a catastrophic creation of quanta of a Kaluza-Klein 
mode in a D-dimensional generalized Kaluza-Klein theory. As the origin of the catastrophic 
creation we have considered a small oscillation of a scale of compactification around a today's 
value of it. To make our arguments definite and for simplicity we have considered two classes 
of models of the compactification: those hj Sd {d = D — 4) and those by Sd^ x Sd2 {di > d2, 
di + d2 = D — 4) , then have shown that the compactification by Sd can be stable for d > 2 
and the compactification by Sd^ x Sd2 can be stable for di > 2. For these stable models 
we have given a hamiltonian for a Kaluza-Klein mode. A 4-dimensional metric and a scale 
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of the compactification have been treated classically, and the Kaluza-Klein mode quantum 
mechanically. The form of the hamiltonian for the Kaluza-Klein mode shows that a so- 
called parametric resonance phenomenon can occur: quanta of the Kaluza-Klein mode can 
be excited catastrophically when frequency and amplitude of the oscillation of the scale 
of the compactification are in a band in their parameter space. We have given formulas 
of a creation rate and a number of created quanta of the Kaluza-Klein mode due to the 
parametric resonance, taking into account the first {p/q = 1) and the second {p/q = 2) 
resonance band. After that by using the formulas we have calculated those quantities for 
each model of the compactification. 

We have shown that for the model of the compactification by 5*^ {d > 2) the first 
resonance cannot occur but the second resonance can occur for the / = 1 mode. For a 
model of the compactification by S^i x 5*^2, if there exists a set of positive integers li and 
I2 satisfying B < 1 then the first resonance can also occur for the corresponding mode 
0Zi«2mim2) where B is given by ( |3.18| ). Finally we have given conditions for the parametric 



resonance to produce sufficient quanta of the Kaluza-Klein mode for each model of the 
compactification. Our conclusion is as follows: (1) in the model of the compactification by 
Sd the parametric resonance does not overproduce the quanta of the Kaluza-Klein mode; 
(2) in the model of the compactification by Sdi x Sd2, if there is a mode satisfying B < 1 
and \B — 1\ = 0{e), then the quanta of the Kaluza-Klein mode can be overproduced by the 
parametric resonance; (3) otherwise in the Sd-^ x Sd^ model, the parametric resonance is so 
mild that the overproduction does not occur. 

As mentioned in Sec. |, if entropy production is negligible, there must not be a catas- 
trophic creation of quanta of a Kaluza-Klein mode. This condition can be regarded as a 
cosmological criterion for compactification. Thus from the results obtained in this paper, 
assuming that entropy production is negligible, we can conclude that the model of the com- 
pactification by Sd-i^ X 5*^2 is ruled out if there is a set of positive integers h and I2 satisfying 
the condition B < 1 and |-B — 1| = 0(e), where B is defined by One may expect 



that infiation followed by reheating can alter the conclusion and save the model. However, 
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maybe it is not the case. In general it is natural that the field o", which represents a scale 
of the compactification, couples to the infiaton field. Hence at the end of the slow rolling 
phase of the infiation we can expect that energy density pa- of coherent oscillation of the a is 
of the same order as energy density of oscillation of the infiaton. In this case the parametric 
resonance of a Kaluza-Klein mode occurs after infiation, if our 4-dimensional universe is 
obtained by the compactification with the compact space Sd^ x and if there is a set of 
positive integers /i and I2 satisfying the condition B < \ and \B — \\ = 0{t). The energy 
density pkk of the created quanta of the Kaluza-Klein mode is expected to be of the same 
order as p^. Hence pKK/Prad = 0(1) is expected after the reheating. Thus, naively, the 
infiation cannot save the model. However, we have to keep in mind that there are, as usual 
in arguments on early universe, loopholes in the above arguments: (1) a late time infiation 
may save the model if there is no enough time after the late time infiation for the ratio 
PkkI Prad to grow Sufficiently; (2) if compactification radius is small compared with experi- 
mental scale but large enough compared with Planck scale then there may be a non-trivial 
range of allowed parameters. To fill up the loopholes is one of the future works we have to 
do. 

Our next question now is whether or not the model of the compactification by Sd and 
the model by Sd^ x S'^j, if there is no integers h and I2 satisfying B <1 and |i? — 1| = 0(e), 
do safely pass the criterion. We can say that the models are not ruled out from our analysis. 
However, unfortunately, we cannot answer this question definitely from only the results 
obtained in this paper since there is a possibility that the structure of the resonance may 
alter seriously when e > 0(1). In this paper we have considered only the case that the scale 
of the compactification oscillates with a very small amplitude (e <^ 1). However, supposing 
that our 4-dimensional universe was realized by a dynamical mechanism of compactification 
from a higher- dimensional spacetime, it is natural that there is an epoch when e ^ 0(1) 
before the epoch we have considered. The former epoch is called broad resonance regime 
and in this regime it is expected that the resonance becomes broad and more efficient. In 
fact, in a simple model of reheating after infiation, it is the case [^]. Thus we have to analyze 
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the former epoch in future works in order to answer the above question. 

Although we have investigated only two classes of models of the compactification, qualita- 
tive properties of parametric resonance of a Kaluza-Klein mode are expected to be common 
in a wider class of models of the compactification. It is because mass of a Kaluza-Klein 
mode and frequency of an oscillation of a scale of the compactification are both expected to 
be of order of the Planck mass in a wider class of models. Hence it is valuable to analyze 
a parametric resonance phenomenon of a Kaluza-Klein mode in more realistic models of 
the compactification with other compact manifolds (eg. Calabi-Yau manifolds) or by other 
ways of stabilization of the compactification (eg. non-trivial configuration of anti-symmetric 
fields, gaugino condensation, etc.). 
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APPENDIX A: FORMULA OF THE FLOQUET INDEX 

In this appendix we consider the following ordinary differential equation for a real func- 
tion Y of z. 



where ^ is a real constant and g{z) is a real periodic function of z with the period 27r. We 
assume that g{z) — g{—z). In this case this equation is called a Hill's equation. For a set of 
linearly independent two solutions Yi{z) and Y2{z), li(27r -|- z) and Y2{2t: -\- z) are hnearly 
independent, too. Hence, by the uniqueness of the solution of the equation, the later two 
solutions are linear combinations of the former two solutions: 




(Al) 




(A2) 



26 



where is a 2 x 2 regular matrix with all components being real constants. Therefore there 
are independent linear combinations Yi{z) and of Yi{z) and 1^2 (-z) such that 

¥,{271 + z) = <f)^Y^{z), 

¥2(271 + z) = <f),Y,{z), (A3) 

where 0i and 02 are eigenvalues of U and are not zero by definition. By introducing a 
constant /i by 0i = e^^'^ , Yi{z) can be expressed by a periodic function (f){z) with the period 
27r as □ 

yi(z) =e^^0(z). (A4) 

Since ( [Al| ) is invariant under the replacement z ^ —2;, Yi{—z) is also a solution. Moreover, 
if 3ft/i 7^ or 2$5yU is not an integer then Yi{—z) is linearly independent of Yi{z). Thus, in 
this case, 02 = e~'^'^^ and the general solution of ( |Al|) can be written as 

Y{z) = Cie'^Xz) + Cse-'^X-z), (A5) 

where ci and C2 are complex constants such that both Y and dY/dz are real aX z = zq |^. 
The constant /i is called a Floquet index. 

The equation ( [Al| ) is said to be stable when is pure imaginary, while unstable when 
/i has a non-zero real part 0. The purpose of this appendix is to obtain a formula of the 
value of /i and to give a criterion when the equation is unstable. For this purpose we use the 
method of averaging [|r^. Let us suppose the case when the function g{z) is parameterized 
by a small parameter e and g{z) = 0(e), and assume that g{z) can be expanded by e. In 
this case by redefinition of the real constant h we can expand g{z) as 

00 00 00 00 

9i^) = E E = 2E E ^'9L'^ cosnz, (A6) 

s=l n=— 00 s=l n=l 



Note that in general the function 0(z) is complex. When /x is real, 0(z) is a real function. 
When /i is real, the coefficients ci and C2 are arbitrary real constants. 
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where (s = 1, 2, ■ ■ ■ ; = 0, ±1, ±2, ■ ■ ■) are real constants satisfying 



gi'^ = 0, 



9':^ = g'-L (A7) 



for all s. 

Let us introduce a parameter A by 



^ \ 2 
P 



h=\^^j + eA, (A8) 

where p and q are mutually prime positive integers. We consider a solution of the equation 
(Al) of the following form 



Y{z) =acos + +Y,e'u^'\a,9,z), (A9) 

\^ / s=l 



where a and 6 are slowly varying real functions of z: 



da °° 



JO oo 

^ = Y^e^B,{a,e), (AlO) 
and u^^^ (s = 1, 2, ■ ■ ■) are functions satisfying 

u^'\a, e, z) = u^'\a, e + 2tt, z) = u^'\a, e,z + 2g7r). (All) 



Substituting the ansatz ( [A9| ) into the Hill's equation (|A1|), we obtain the following 0(e) 
and 0(e^)-equations. 



oo ^2 „2 



E ^^u^:\a,e)e-/^ = a{ £ + ^i?i(a, ^^) - A ) cos + ^ 

i=-oo ? \n=-oo 9 / \Q J 



Note that in order to obtain the formula of the value of we only have to seek a special solution 
of (Al) because of the form of the general solution (A5). So we concentrate on this form of the 
solution. 
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+—Ai{a,9) sin (-z + 9 
(1 \(1 ) 



(A12) 



l= — oo ^ 



+ {^JLa, + a^^, + .^B, + 2^,B,j sin + flj 

where the fourier components M^^^(a,6') and v!)^{a,Q) of ^'^^^(a, 6*, 2;) and u^'^\a,9,z) have 
been introduced: 



inz/q 



(A14) 



and ^(2) is defined by 



[ f; ^«e-^-A]w«(a,^,z) + a f] ^(V-cos(^^ + ^ 

\n=— 00 / n=—oo \^ ) 



(A15) 



From the 0(e)-equation ( |A12|) , we obtain the following form of A\, B\ and {n 7^ 



1 



2p/g 



A cos 2^ 



(A16) 



and 



f , (1) 



-id 



(A17) 



Without loss of generality we can set to be zero by redefinition of a and 9. Thus 



u 



'"(0,9, 



2)2 — ^2 



(AI8) 



18 When n is not an integer we define (7!*^ by Qn^ = 
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Next, from the O(e^) -equation (|A13|) , we obtain 



2p/q \ da 



') a ■ 2p/q 



e - g-LpC 

, -(2) ^ 



2p/q 



where we have introduced the fourier component g^:"^^ of g^'^^ by 
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Substituting ( |A16| ) and ( |A18| ) into (|A19|) , we obtain 
A2(a,e) 



2p/q 



■e^ yn y2p/q-n (2) 
^^O^/^ IT ^ 32p/q 



n=i0,2p/q 



n{2p/q — n) 



sin 2^, 



g(l)2 



i2p/q) 
1 

2p/g 



yn y2p/q-n (2) 

^ n{2p/q - n) ^^p/g 



cos 2^. 



_n7^0,2p/ij 

Finally we seek the formula of the value of fi up to order O(e^) by using the 
obtained here. Let us introduce slowly varying functions x and y of z hj 

x{z) = a cos 9, 
y[z) = asm 9. 

The following equation for x and y can be obtained from (|A10|) , ([A16| ) and (|A21|) . 



d 

dz 



where \^ is a 2 x 2 matrix of the form 











\ 


x{z] 




x{z] 






= V 








) 






J 



V 



2p/q 







-(<5 + A)^ 



V 



-{5 -A) 



The constants 6 and A are defined by 
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A = A + e 



\ - y-n y2p/q-n (2) 



(2p/g) 



g(l)2 



2\ _ y 



(A25) 



Since up to order 0{e'^) the Floquet index yU can be understood as an eigenvalue of V , it is 
given by 



/i = /i± 



52- A2 + Ore3). 



(A26) 



Thus the necessary and sufficient condition for the equation (|A1|) to be unstable is 



-|<5|<A<|5|. 



(A27) 



This inequality can be solved with respect to A by using the fact that A = 0(1). The result 
is as follows up to 0(e). 



Ao-|(5|<A<Ao + |(5|, 



(A28) 



where Aq is defined by 



9. 



(1)2 



(A29) 



APPENDIX B: NUMBER OF QUANTA EXCITED BY PARAMETRIC 

RESONANCE 

In this appendix we derive a formula of number of quanta excited by the parametric 
resonance. Although the formula can be found in literature 0, we derive it in this appendix 
for completeness. 

For example, see Ref. Q. However, their derivation is incomplete as will be stated later. The 
obtained formula of the number of quanta is exactly same in Ref. [^] and in this appendix. 
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We consider a quantum mechanics of an oscillator described by the time dependent 
Hamiltonian 



1 r 



H =- 
2 



p2 + n\t)Q' 



(Bl) 



where VL{t) is a positive function of time t, Q and P are operators representing a coordinate 
of the oscillator and its conjugate momentum: 



First let us define the time dependent operator a{t) by 

J f ndt 



(B2) 



a{t) 



nQ + iP) . 



(B3) 



The operator d{t) and its hermitian conjugate have the standard commutation relation 
of creation and annihilation operators: 



a, a 



1. 



(B4) 



Since the hamiltonian can be expressed as 



2) ' 



(B5) 



the time dependent vacuum state |0t), which is defined by a|0() = 0, minimizes the value 
of the Hamiltonian at a moment t. Hence, the creation and the annihilation operators d{t) 
and a){t) regard, so to speak, instantaneous quanta at the moment t. 

The equation of motion for Q and P is equivalent to the following differential equation 
for the operator a{t) Q: 



"-2n" 



A general solution of this equation is written as 



(B6) 
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Hereafter the dot denotes the derivation with respect to the time coordinate t. 
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a{t) = a{t)a{0) + f3*{t)a^ (0), 
where a(t) and (3{t) are complex functions satisfying the differential equation 



(B7) 



(B8) 



and the initial condition 



a(0) = 1, /5(0) = 0. 



(B9) 



The number of the instantaneous quanta at the moment t for the vacuum state |Oo), which 
is annihilated by a(0), can be expressed by j3(t) as 



iV(t) = (Oo|at(t)a(t)|Oo) = |/3(t)r 



(BIO) 



By using (p3| ) and (pTl), the operator Q is expressed by the operator a(0) and its her- 
mitian conjugate a"''(0) as 



g(t) = gW(t)a(o) + g(-)(t)at(o) 

where Q^^\t) and Q^^\t) are defined by 



(Bll) 



or equivalently, 





1 




Q^-\t) = 


1 








a{t) = 






-i f ndt 


m = 





Qdt 



(B12) 



r](t)QW(t)+zQW(t)" 

'n{t)Q^-\t)-tQ^-\t) 



(B13) 



Since a(0) and a^(0) do not commute with each other, the equation of motion of Q{t) means 

+ fi2Q{±) = 0. (B14) 
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It is shown that the initial condition ( [B9| ) is equivalent to 

1 



Q(-)(0) 
by using the relation 



2fi(0) 



(-)i 



(B15) 



i J ndt 

Pit) = ^j== {n{t)Q^-\t) + tQ^-\t)) . 



(B16) 

Thus, in order to calculate the number of instantaneous quanta N{t), what we have to do 
is to seek the solution Q^^\t) of the differential equation (P14|) with the initial condition 
( P15|) , and to substitute it into (P10|) with (|B13| ). Note that the function Q^~^ is complex 
while the classical coordinate Q corresponding to the operator Q is real. 
Let us consider the case when the function Q{t) is given by 



n^{t) =uj^[h- g{iut)], 



(B17) 



where u and h are positive constants, g{z) is a real function given by ([A^). In this case the 
corresponding differential equation for Q^~^ is the Hill's equation with Y replaced by 
Q*-^-*, provided that 



z = uit. 



(B18) 



In Appendix ^ we have obtained a real solution of the form (|A9|) . We write the real part 
and the imaginary part of Q^~^ in the form (1^4 9|) , respectively: 



P 



P 



Q( >{t) = {xi + 1x2) cos yj^^j ~ + *?/2) sin y-^j + 0{e), 
where Xi^2 and yi^2 are real functions satisfying the following differential equations. 



(B19) 



dz 





\ 






\ 


Xl,2{z) 




Xl,2{z) 




= V 






) 






J 



+ 0{e' 



(B20) 



where the matrix V is given by (|A24|) . For the functions xi^2 and yi^2, the following initial 



condition in the lowest order of e is derived from the initial condition ( |B15| ): 
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^i(0) = 2/2(0) 



X2(0) = 2/i(0)=0. 

Hence, in the lowest order, the solution of ( |B20|) is given by 

1 



Xi 



yi 



X2 



y2 



uj ■ 2p/q 
1 



: cosh n+z, 



UJ 



■ 2p/q\5-A 



^ smh n+z, 



Since it is shown from ( P13|) and (|B19|) that 



1 5-A . , 

\ ^sinh/i+^, 

. 2p/q + 
: cosh jjij^z. 



m)\ 



UJ ■ 2p/q 



(xi- 2/2)2 + (a;2 + 2/i)2 + 0(e) 



the number of quanta N{t) defined by (piOj ) is |^ 

1 



N{t) 



2 sinh^ jjL+z + 0(e). 



(B21) 



(B22) 



(B23) 



(B24) 



1 - (A/5J 

Finally we consider the case when the value of 5 and A are not constants in time but 
they are given as functions of z. We assume that their change is so slow that 

d 



— ln|5± A| 

dz 



< (B25) 
In this case /i+z in (|B22| ) is replaced by / fi+dz. Thus the number of quanta in this case is 



N{t) ^ sinh^ (^J fx+dz^ . 



(B26) 



^-^ In Ref. they used the equation /3 = + 0(e) to derive this formula of N{t). However, 
the term //^/3 is of order 0{e'^) since fj,^ is of order 0(e). Thus their derivation is incomplete. 
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